HTC-AD 71-59 


HUGHES 7. 62MM EXTERNALLY Е 
ARMOR MACHIN 


PROGRESS REPORT 3 
8 March 1971 Through 9 June 1971 


Cont >t DAAF03-70-C-0011 


10 June 1971 


Digitized by: 


SMALL ДЕМІ)! 


MALLARMSOFTHEWORLD.COM 


HTC-AD 71-59 


HUGHES 7. 62MM EXTERNALLY POWERED 
ARMOR MACHINE GUN 


PROGRESS REPORT 3 
8 March 1971 Through 9 June 1971 


Contract DAAF03-70-C-0011 
10 June 1971 


әзәтавт}цолд 


НТС-АР 71-59 


HTC-AD 71-59 


TABLE OF CONTENTS 


Page 
INTRODUCTION 5-40-00 ка Са ОЛОВО adi ov Fea И РОДЕ Y ЛЕКТ. 1 
ЭУММАНҰ е бы sheen eye EEEIEE е ее Lene T 3 
SYSTEMODESIGN/O оа ен аа Қы, евра з DOCE EE 5 
COST REDUCTION ...... ілуде ы» Wate Ша ue 8245646 9 
МАФНАМӘРЕІ ИУ 22072 e а еа басенот $30 gia tr тойлар 6 ова ее T t 10 
АРРЕМОГ ве pok deia ЕЕ qi ава дода, т в e s 12 

iii 


HTC-AD 71-59 


INTRODUCTION 


The Externally Powered Armor Machine Gun (EPAM) shown in the frontis- 
piece and the proposed layout (Figure 1) is being designed and developed and 
a model is being fabricated under Contract DAAF03-70-C-0011 with U. S. 
Army Weapons Command. The primary characteristics of the EPAM design 
are: 


Externally powered 

Dual feed 

600 shots per minute 

Cam operated 

Single barrel, quick replacement 
7.62mm NATO ammunition and M13 link 
Compact size, the same as M73 
Sprocket feed 


While this initial design is scaled for 7.62mm NATO, the general concept is 
adaptable to a range of calibers through 30mm. 
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Figure l. Layout - Externally 
Powered Armor Machine Gun 
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SUMMARY 


The Third Contract Review meeting was held at Hughes Tool Company - Aircraft 
Division, Culver City, California, on 8 March 1971. The second progress 
report was delivered and a review was made of work accomplished to that date, 
Specific items discussed included the following: 


Status of system design 

Level of stress analysis required on model 

Necessity for spring packs on bolt carrier 

Extent of breadboard testing 

Complete review of system math model 

Changes and additions to system math model 

Possibility of a 10-percent cost overrun resulting from a change in 
HTC-AD's general engineering overhead 


For the 8 June 1971 Contract Review, the System math model has been com- 
pleted and now includes all changes requested during the previous review. Тһе 
gun system layout is complete, as are all details of the bolt, motor/cam drive, 
and significant portions of the feed systems. The bolt has been released for 
fabrication and the cam has been completed (Figure 2). 


In addition, an investigation has been made of the contractual obligations and 
remaining funds in an effort to reduce the impact of the increase in the general 
engineering overhead and G&A rate. The object of this study is to identify the 
most significant portions of the scope of work that can be completed within the 
funds available. 


During the coming period, it is anticipated that the following work will be 
accomplished: 


e. Complete detail design 
e. Initiate fabrication of all components 
e Assemble bolt, cam, and drive assembly 


Figure 2. 


Bolt Drive Cam 
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SYSTEM DESIGN 


System design is proceeding extremely smoothly and on schedule. The final 
layouts have been completed and detail design is 40 percent complete. The cam 
and bolt group have been released for fabrication. 


А number of representative drawings are included in this report (Figures 3, 4, 
and 5) to illustrate some of the design work accomplished during this period. 


Of particular interest is Figure 3, which shows the proposed conformation of 
the model's cam and drive system. А hand crank system of operation is 
shown; the installation of an electric motor has been eliminated in order to 
conserve funds but is shown in phantom line. Also shown in this drawing is 
the simplified construction by means of flat plates and cap screws that is being 
used to reduce construction costs. 


Included in the Appendix to this report is a stress analysis of the firing pin 
spring and bolt that have resulted from the redesigned bolt system. This 
system allows for significantly lower firing pin stress, easier disassembly, 
and stronger construction than that originally proposed. 
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COST REDUCTION 


An increase in the Company Engineering overhead and G&A rates effectively 
reduces the remaining funds available for completion of the design and fabrica- 
tion of the gun model on this contract. To absorb the effect of this reduction -- 
thus avoiding а 10-percent contract overrun -- and still permit the completion 


of the important aspects of the Scope of Work, certain changes in the program 
are recommended, as follows: 


1. Forward ejection will be eliminated from the model and in its 
place a simple side-ejection system will be utilized. This is not 
necessarily a permanent change, and forward ejection can be 
reinstated at a later date if funds were to be made available. 


2. The model will be hand powered by a crank rather than electrically 
driven. Provisions can be made for the later addition of an electric 
motor and control box. 


3. АП planned travel to WECOM by HTC-AD personnel on this contract 
will be eliminated and those travel funds will be applied to the 
engineering task. 


4. The overall construction of the gun model will be simplified and 
such niceties as simple "takedown"! have been eliminated. The 
essential elements of the EPAM concept will be clearly illustrated, 
but nonessential items will not. 


5. Scheduled breadboard testing will be eliminated and all testing will 
be done at the completion of the model. 


6. Reports will be simplified. 


7. It is recommended that one review meeting and report be eliminated 
and that the program schedule be revised as shown in Figure 6. 


With these changes and an anticipated reduction in the estimated cost of model 
fabrication, it is hoped that the program can be completed within the available 
funds, 
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MATH MODEL 


The mathematical model of the EPAM system has been completed and the 
changes and additions requested during the last program review have been 
incorporated. In view of the funding situation, no further expansion of the 
model is planned. 


А complete and corrected version of the math model is included in the Appendix 
to this report. 
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HUGHES TOOL COMPANY-AIRCRAFT DIVISION 
macys Luring Fon Energy MODEL REPORT NO. PAGE 
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Spring 365 OD 13 Total coils ends closed 4 ground // асўуўе coils 
-,055 WD, solid М, 770 before grind ,720 After grind 

1310 М.р, rales 40,5 у 
AF: WH,-WHa-.406 мире Ба E) + 40 (АР) = 158% 
ДР = 406A) х 105 И (vate) = 16,42" | Free H 2 В + WH 
vate 
Р, = 13.8 (оуР)-/0,42Ҡ% АР) = 5,59% dup +, 7815 |, 825 
2 40,5%» 
Py = 13.9*Cau.P) + №.42(АР) = 22,01 
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y= 22.01 


$4= /33,200 RE а ҚС 
аз 
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$ = .450 
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55“ 112,000 
МНз.875- 
Ignition 
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CHECKED ву. 


force = ,225 v* fT у 80,000 Ps, = 12, 700% 


Sheay ауса = ‚4906 X, 82x 3 =, 389 57 in 
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material (Carpenter 58) shear strength = 112,000 pst 


Safety factors 8441 ( proof rounds ) 
505%) Cstandard Ammo) 


Comments | lug will hold max proof round 
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EPAM MATH MODEL 


A computer math model to simulate the Hughes 7. 62mm Externally Powered 
Armor Machine Gun (EPAM) has been formulated and programmed as part of 
this effort. The math model utilizes the Tymshare, Inc., Super Fortran 
language (IBM LevelH Fortran IV). 


DEFINITION OF SYMBOLS 


a = acceleration of belt 

b - moment of belt 

c - moment arm for rotating track reactions 

d = distance, cg of bolt to center of cam roller surface 

e 7 distance from roller center to nearest track contact point 
Cr = coefficient of restitution for bolt-carrier impact 

F = driving force 


Fa = axial inertial force of bolt and round or of bolt and case 
E 7 force required to drive ammunition belt 
lequiv = mass moment of inertia of all rotating parts 
= length of bolt travel 
M = МВ + МС = mass of bolt unit + mass of carrier 
MB = mass of bolt 
Мс = mass of carrier 
N = normal force on roller 
Ма = axial component of the normal force of the roller 
№ = transverse component of the normal force of the roller 
R - radius of feed sprocket 
Re = cam radius 
В = frictional resistance due to track reactions 
Ry = track reactions due to rotational forces 
Rt 7 track reactions due to tipping forces 
Tcam = torque about gun axis 
VB = axial velocity of bolt 
Ус = axial velocity of carrier 


WTbelt = total weight of ammunition belt 
У = axial displacement of bolt 

Y - axial velocity of bolt 

Y = axial acceleration of bolt 

a = angle of cam path (slope) 

8 = angular displacement of rotor 
8 = angular velocity of rotor 

8 - angular acceleration of rotor 
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u = load factor for ammunition belt 
Ur = coefficient of rolling friction 
Ut = coefficient of friction of track 


MODEL OPERATION 
Basically, the model works with the following equation of motion. 


EQUATION OF MOTION 


(ear + Нееа + редеа) 9 = Tmotor - Team - Tfrict - Тъей (IA) 
where 

Icam = moment of inertia of cam 

Igear = E moments of inertia of gears weighted for any reduction 
factor 

Iteed 7 moment of inertia of feed sprockets and rounds adjusted 
for any reduction factors 

5 = angular acceleration of cam 

Tmotor 7 motor torque as a function of rpm 

сака 7 torque due to the loads on the cam, including friction and 
bolt inertia loads as a function of the cam curve 

Trrict = torque due to gear friction 

Thelt = torque due to weight of rounds in belt being moved 


Solving for the angular acceleration, 8, yields 


T -T -T -T 
- motor cam frict belt 
8 = (IB) 
I EA +1 
cam feed gears 


NUMERICAL INTEGRATION METHOD 


The integration technique employed by the program is Adams four-point 
method. The integration formula is as follows: 


уч (55Y,-59 Y, +37 У а ан (ш) 


In order to obtain nts а Y ge and Ya. to start the integration 


$ 
process, a Taylor series expansion is performed as follows: 
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M yg *àty, 0 ) 
EP TP + у) e) 

a = vote (2 о. n Ы зу) ey 

yi) E ye on _ ки (4) (ш) 
9 z ,e £x (oy а Е + уй (i Қ s, Фф y) (7) J 


Equations (3), (5), and (7) are then used as solutions for the first three 
points. 


The math model uses the following sequence in integrating: 


= ir (Equation of motion, IB) 


ntl 
TIME = TIME + ATIME 
ntl n 
Ber 17 8, + | 8, dt (Equation II above) 


Gna. = ба + ее dt (Equation II above) 


MOTOR CHARACTERISTIC CURVE 


The math model utilizes a torque versus rpm curve to simulate the drive motor. 
Upon attaining steady-state rpm, the motor is "shut off' and the program 
simply calculates the torque that would be required of the motor to keep the 
system running at steady-state rpm. Ап exception is made at the point where 
the carrier picks up the bolt. The resulting velocity after impact is decreased 
substantially so that the motor must operate at full power in order to restore 
the system to steady-state rpm. 
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CAM CURVE 


The cam drive used in the math model is a modified-trapezoidal type, with 
41.25 degree constant velocity portions, The constant acceleration portion 
(parabolic) has the following characteristics 


2 


cam 


where A = 1.6609/ft = the curve constant. 
cam 


The cycloidal portion of the curve is defined by а quarter wave with the 


peak point equal to Асау. А layout of the curve is shown on page 21 of this 
section. 


CAM FORCES 


The torque on the cam (TCAM) is obtained as follows: 


iD = NR 
cam CRC 
where 
м, = the transverse component of the cam normal force 
R = the cam radius 
с 
Also, 
м, = Мапа tu, №05 а 
where 


N and « are as defined in the free-body diagram of page 23, and 
uy N is the resistance induced by rolling friction 


Ма = М соѕ а -ur N sing 
а r 


The track reactions due to rotational forces are found by balancing moments 
in the plane perpendicular to the bolt-carrier axis. 


eR, = d[F - F,) 


ва (2) (ғ - ғ.) 
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ACCELERATION CURVE 
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Figure 7. Cam Characteristics 
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The track reactions due to tipping forces are found by balancing moments in 
the vertical plane parallel to the bolt-carrier axis. 


The frictional resistance due to track reactions is 
Rey = +2, (В, + қ.) 


Ry, has the same algebraic sign as Y. 


The normal force on the cam roller is found by balancing the axial forces, 
thus У Еу= 0 


Ка +, М-М +R, = 0 


fr 


Substituting all the expressions containing N in the above equation yields 
e d d Я 
N= |е. Los +. (1+) +2 нь - 2% ugue | sin a 


e d d 
“рын (15$) зін дү ee бре то Јана 


The corresponding equation for the deceleration phase becomes 


d 
М= [Fa || ем (не) +22 an, - 224, - 1.0 сов o 
е d d А 
[ut (1+=) нЕ ue] sin а] 
Letting 
e d d 
€ rupta (1+5) +2 ы-2 нуш 
~ е а 4 
Сазы, (1+2) +2 ПЕН др мет 1 
е 4 
€, = -fur tut (1+2) + 22 це + 25 urut 
é d 
бутика (14$) +25 шш, 25-1 
Ка = МУ +K, У, + Kos Ура + Fetrip 
22 
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yields 
N= Е - 5 / ша 
N=- (мук, аға жы БІ (с sina + C, соз a) 
(The negative sign is added here to cancel the sign of С», which is negative) 


For the deceleration portions of the cam curve, 


Ме Еу" - (МУ - F strip )/ 6 sin a+ C, cos о) 


where 
M = mass of bolt-carrier combination 
Y = axial acceleration of M 
5 - combined spring rate of carrier spring pack and cartridge 
(used only when simulating the jamming of a round in the 
chamber) 
Y. = the deflection associated with Ks 
K = spri te of firi i 
29 pring rate of firing pin 
Yom = the deflection of the firing pin spring 
Е ое the average load required to strip а round from the link 
strip 


= 6.8 Ib (test data) 
It should be noted that the terms comprising F, above act only during those 


portions of the cam curve as defined by the corresponding input cam positions 
(see input data section). 
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BOLT-CARRIER IMPACT LOADS 


At the indicated cam curve position, THIMP & THREL (input data, items 21 
and 22), the carrier releases the bolt (THREL) and picks it up again during 
its rearward motion (THIMP), at which time a momentum exchange takes 
place between the moving carrier and the stationary bolt. 


Three simplifying assumptions will be made at this point. First, the bolt 
and carrier will be assumed to be moving at the same velocity after the 
initial impact, thus ignoring any subsequent minor collisions that take place 
before the two are locked together. Second, the velocity of the system (bolt, 
carrier, plus cam and gears) after impact will be given by that of the first 
mass (carrier, cam, and gears) after impact. Finally, the collision is 
assumed to take place during the rearward acceleration portion of the para- 
bolic curve. 


The momentum of the first mass (carrier, cam, and gears) before impact is 
given by 


S 6 sina соба 
"egy 200 ов 
Momentum = Ме Ус + 


Вс 
where 
Ме - mass of the carrier 
Re = radius of cam 
Ме = axial velocity of the carrier 
lequiv = combined moment of inertia of the cam and gears 
а. = angle of cam slope 
8 - angular velocity of the cam 


For the parabolic portion of the curve, the carrier velocity and 8 are 
related as follows: 


dv 2 „аде 
Verna = 2А 


= the parabolic curve constant 
the angular displacement along the parabolic portion of the 
curve 


>» 
o 
ы 
3 

"a 


Substituting yields 
+ equiv Ve 


2Acam Вс 


1 ,.,Sinacosq 
BC Ст от, 


Momentum - Ме 
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' 
The velocity of the system after impact, Ve, is then 


sina cos а, 


equiv 
Mg + S889 | Ve + С, Mp Ус - Vp) 
қ БА сата Ке 8 
Vere I, wn sina cos о. 
Me + 3 * Mp 
2Acam Rc 8 
where 

e 7 the coefficient of restitution for the impacting masses 
Ма = the mass of the bolt 
Ув = the velocity of the bolt before impact (v4 = o) 


BELT TORQUE 


R 
NT кш. (RADSPR) 
ж 
ED. ut 5 


From the free-body diagram, the force, F., is given as 


1 


belt 
By PH WT git a © 
RE 
Fit WT еј (ya xi) 
where 
u = 1.8 = load factor for an ammunition belt in a horizontal 
chute (Reference 2) 
WT. = weight of the ammunition belt 


belt 
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2 
32.17 ft per вес - gravitational constant 


а 
" 


a zm аб the acceleration of the belt 


ә: 
" 


the angular acceleration of the cam (since the feed sprocket 
rotates at one-twelfth the speed of the cam, 8/12 is the 
angular acceleration of the feed sprocket) 


Hence, the torque due to belt loads is 


RE 
Traine Soa PT ој (: неа * RE) 
where 
Ge aia = the G-loading under which the belt will be moving. 
27 
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The following is a Fortran name list and description of all the inputs in the 


order they should read. 


Data 


Location 


1 


10 


11 


14 


15 


Fortran 
Name 


ACAM 
GFIELD 
RADIUS 
RADSPR 


CR 


OMEGA 


MB 
WTBELT 
LCAM 


IGEAR 


IFEED 


ICAM 


KP 


КТ 


Description 


Cam curve constant 
Gravitational field 
Mean cam radius 
Radius of feed sprocket 


Coefficient of restitution for bolt- 
carrier impact 


Steady state angular velocity of 
cam 


Mass of carrier 

Mass of bolt 

Weight of ammunition belt 

Axial length of cam 

Moments of inertia of all gears 
normalized to the angular velocity 
of cam 

Moment of inertia of feed sprocket(s) 
normalized to the angular velocity of 
cam 


Moment of inertia of cam 


Spring rate of carrier spring rack 


Damping coefficient for gears 


28 


Units. 
1/Ft 
G's 

Ft 

Ft 


None 


Rad/Sec 


Slugs 

Slugs 

Lb 

Ft 

Ft- Lb-Sec? 
Ft-Lb-Sec? 
Ft-Lb-Sec? 
Lb/Ft 


Lb-Ft 
Rad/Sec 
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Data Fortran 
Location Name Description Units 
16 KR Spring rate of rounds Lb/Ft 
1T KPIN Spring rate of firing pin mechanism Lb/Ft 
18 YINIT Bolt position when jamming begins Ft 
19 YPRE Preload of carrier spring pack 
20 ASIN Length of cycloidal portions of cam Rad 
(NOTE: Allangular cam positions 
are referenced with respect to the 
rearmost follower position) 
21 THIMP Angular position of cam when Rad 
carrier impacts bolt 
22 THREL Angular position of cam when Rad 
carrier releases bolt 
23 THTCRU Angular position of cam when Rad 
jamming begins 
24 THRAM(1) Angular position of cam when Rad 
cartridge ramming begins 
25 THRAM(2) Angular position of cam corre- Rad 
sponding to end of peak stripping 
load 
26 THRAM(3) Angular position of cam when Rad 
cartridge ramming terminates 
27 THPIN(1) Angular position of cam when Rad 
contact is first made with firing pin 
Spring during forward motion of bolt 
28 THPIN(2) Angular position of cam during for- Rad 


ward motion where firing pin 
becomes fully cocked 


29 


Data 


Location 


29 


30 


ЗА 


32 


33 


34 


35 


36 


31 


38 


39 


40 


41 


42 


Fortran 


Name 


THPIN(3) 


TH(1) 


TH(2) 


TH(3) 


TH(4) 


TH(5) 
TH(6) 
TH(7) 
FSTRPI 


FSTRP2 


FLOCKI 


UR 


UT 


Description 


Angular position of cam where bolt 
locking terminates and firing pin is 
released 


NOT USED 


Angular position of cam when first 
parabolic portion of curve termi- 
nates 


Angular position of cam when first 
cycloidal portion of curve termi- 
nates 


Angular position of cam where con- 
stant velocity portion of curve 
terminates 


Angular position of cam where 
second parabolic portion of curve 


terminates 


Same as TH(1) + 360% 


Set 


Same as TH(2) * 360 танана у 


Same as TH(3) + 360° 
Initial peak stripping load 


Average stripping load excluding 
initial peak load 


Average load during bolt unlocking 
phase 


Coefficient of rolling friction for 
cam follower 


Coefficient of sliding friction for 
carrier track 
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Units 


Rad 


Rad 


Rad 


Rad 


Rad 


Rad 
Rad 
Rad 
Lb 


Lb 


Lb 


None 


None 


Data 


Location 


43 


44 


45 


46 


47 


66 


67 


86 


Fortran 
Name 
BB 
CB 


CB 


EB 


RPMT(1) 


ВРМТ (20) 


ТТАВ(1) 


TTAB(20) 
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Description Units 


Moment arm of tipping track In. 
reactions of carrier 


Moment arm of rotational track In. 
reactions 
Distance from bolt-carrier center In. 


of gravity to cam roller surface 


Axial distance of carrier between In. 
rearmost edge of track and contact 
point on roller surface 


First value of rpm for the motor RPM 
rpm versus torque table in ascend- 
ing rpm order 


Last value of rpm for table RPM 


First value of torque for rpm versus Oz-In. 
torque table 


Last value of torque for rpm versus Oz-In. 
torque table 


The following free-format input data precedes data items 1 through 84: 


None 


None 


None 


None 
None 
None 


None 


ND 


DT 
x 


рх 

TFINAL 
KPRINT 
CYCLES 


Total number of data items None 
(1 through 84) to be read in if 
different from built-in block data 


values 
Integration step size Sec 
Initial angular position of cam Rad 


(0 is rearmost dead center) 


Initial angular velocity of cam Rad 
Final program time Sec 
Print interval None 
Number of firing cycles to be run None 
after attaining steady-state RPM 
(OMEGA) 
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Fortran 


Name 


TIME 


THETA 


DTHETA 


DDTHETA 


TORQUE 


POWER 


RPM 


FN 


TBELT 


DDY 


ALPHA 


THTA 


OUTPUT DATA 


Description 


Counter, for internal program bookkeeping only 
Time from application of power to the cam 
Total angular displacement of cam 

Cam angular velocity 

Cam angular acceleration 


Motor torque available during start-up, re- 
quired during steady state 


Motor power available during start-up, re- 
quired during steady state 


Cam angular velocity 

Normal force between cam and cam follower 
Torque on system due to ammunition belt loads 
Axial acceleration of the bolt-carrier system 
Time dependent value of the cam curve slope 


Angular distance measured from the origin of 
the parabolic portions of the cam curve 
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Rad/Sec 
Rad/Sec? 


Lb-Ft 


HP 


RPM 
Lb 
Lb-Ft 
Ft/ Sec? 
Rad 


Rad 


EPAM MATH MODEL OUTPUT 


>RUN 
Os 6001+ +0, 60s 64952 


INITIAL CONDITIONS 


TIME STEP SIZE FINAL TIME 
0.000 +0010 +4000 
MOMENTS OF INERTIA CAM BOLT 
CAM GEARS FEED SYS RADIUS MASS 
FT-LB-SEC2 FT SLUGS 
•0008700 «0000277 «0000017 411717 „09969 
I. TIME THETA — DTHETA DDTHETA TORQUE POWER 
SEC DEG ^ RAD/SEC RAD/SEC2 LB-FT HP 
1 0.000 0.00 0.00 0.0 5.833 0.000 
2 «004 2.90 25.18 6153.4 5.666 ‚259 
з  .009 14.35 54.11 — 5360.1 5.082 “500 
а «014 33.42 78.04 — 4127.8 4.317 +612 
5 +019 58,39 95.07 2696.8 3.623 4626 
6 «024 87.24 105.16 1359.6 3.139 +600 
7 +029 118.07 109.38 579.8 2.878 +572 
8 «034 150.19 116.12 1629.9 2.520 +532 
9 .039 184.55 123.41 1374.5 2.060 +462 
10 •044 220.23 127.33 810.7 1.810 +419 
11 .049 258.11 131.55 2474.9 1.510 +361 
12 .054 297.69 143.02 1187.9 +702 4183 
13 .059 339.40 146.79 381.9 “405 “108 
là 4064 381.55 146.94 -307.7 +395 +105 
15 +069 423.30 144.13  -769.6 “614 +161 
16.076 469.72 101.90 2265.4 3.342 +619 
17 +081 499.78 108.50 1942.8 2.926 4577 
18  .086 532.20 117.37 1643.2 2.445 +522 
19 .091 566.91 124.71 1293.1 1.978 +448 
20 .096 603.77 134.24 2489.2 1.297 4317 
21 2101 643.87 144.85 1666.2 +556 +147 
22  .106 686.30 150.43 546.1 +163 +045 
23 .111 9.51 150.29 -578.1 +172 +047 
24  .116 51.94 145.27 -1326.7 “5923 “138 
25 „121 92.53 137.92 -1567.8 1.053 *264 
26 +126 130.98 131.69 -87.4 1.499 • 359 
27  .131 169.21 135.34 738.7 1.217 “300 
28 .136 208.30 136.34 117.8 1.154 +286 
99 .141 247.87 142.74 1756.3 *725 4188 
30 .146 289.79 149.05 846.3 4255 +069 
31  .15! 332.93 151.45 143.4 +095 +026 
32 .156 | 376.26 150.52 -501.2 “158 “043 
93 .161 | 418.90 146.78 -955.2 +405 “108 
34 ч168./ 466.23 104.61 1931.7 74 “604 
35 .173 496.89 109.89 1652.3 2.851 +570 
36 .178 529.74 118.31 1632.0 2.385 513 
37 .183 564.66 125.68 1237.6 1.916 „438 
38 .188 601.78 134.13 2479.3 “306 “318 
39 .193 641.85 145.08 1679.3 .538 “149 
40 .198 684.36 150.76 568.9 “142 +039 
41 «203 7.68 150.76 -562.6 4141 +039 
42 — .208 50.25 145.76 -1333.1 m ‚128 
43 .213 90.98 138.37 -1575.7 1.025 «258 
44 +218 129.50 131.84 -261.7 1.488 +357 
45.223 167.69 135.54 685.7 1.205 .297 
45 „227 199.02 135.86 360.4 1.184 4292 
AVERAGE TRANSIENT POWER = .54639996 НР 


AVERAGE STEADY STATE РОМЕН = 
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FRINT INTERVAL 


5 


RPM 


0.0 
240.4 
516.7 
745.2 
907.8 
1004.2 
1044.5 
1108.8 
1178.5 
1215.9 
1256.2 
1365.7 
1401.8 
1403.2 
1376.3 
973.1 
1036.1 
1120.8 
1190.8 
1281.9 
1383.2 
1436.5 
1435.2 
1387.2 
1317.0 
1257.6 
1292.4 
1301.9 
1363.0 
1423.3 
1446.3 
1437.3 
1401.7 
998.9 
1049.4 
1129.7 
1200.1 
1280.8 
1385.4 
1439.6 
1439.7 
1391.9 
1321.3 
1259.0 
1294.3 
1297.4 


BOLT 
IN 
0.000 
“001 
+017 
“093 
+284 
+634 
1.162 
1.837 
2.577 
3.329 
4.092 
4.633 
4.922 
4.918 
4.623 
3.954 
3.328 
2.646 
1.899 
1.126 
“483 
“095 
“008 
+225 
nm 
1.430 
2.246 
3.088 
3.909 
4.546 
4.896 
4.935 
4.668 
4.016 
3.395 
2.699 
1.948 
1.165 
+509 
+106 
+005 
„210 
«690 
1.398 
2.214 
2.888 


TIME 

SEC 

0.000 
•004 
+009 
+014 
“019 
„094 
+029 
“034 
+039 
“044 
“049 
“054 
+059 
“064 
“069 
+076 
“081 
+086 
4091 
“096 
“101 
“106 
elll 
4116 
„121 
•126 
4131 
.136 
4141 
146 
“151 
+156 
+161 
“168 
4173 
“178 
+183 
4188 
+193 
+19в 
+203 
«208 
“213 
„218 
„227 


FN 
LB 
«000Е%01 
4135Е%01 
4629Е%01 
*134E*02 
«205Е%02 
.264Е+02 
.278Е+02 
• 105Е+02 
*765Е+01 
+ 108E+02 
+ 149E+02 
+148Е+02 
• 166Е+02 
•19гЕ+02 
.208Е+02 
*152E«02 
+ 122E+02 
+9 42Е+01 
•752Е+01 
•173Е+02 
«245Е%02 
„28 6Е+02 
«З324Е%02 
“330Е%02 
«331Е%02 
•171Е+02 
«487Е%01 
• 103E+02 
+165Е+02 
•165Е+02 
+177Е+02 
+198Е+02 
.211Е+02 
.172Е+02 
4144Е%02 
+882E+01 
“737Е%01 
+ 169E+02 
*245Е+02 
+28 7E+02 
• 325Е+02 
*331E*02 
«332Е%02 
+ 189E+02 
«807Е%01 


TBELT 

LB-FT 
%000Е%01 
«937Е-01 
«899Е-01 
«838Е-01 
“761Е-01 
*687E-01 
*630E-01 
+702Е-01 
+675Е-01 
%629Е-01 
«590Е-01- 
*675E-01 
*689E-01- 
*590Е-01 
%562Е-01 
•502Е-01 
+69 0Е-01- 
«698E-01 
«613Е-01- 
*731E-01 
*701E-01- 
*642E-01 
%580Е-01 
+533E-01 
*517E-01 
«555Е-01 
“642Е-01 
«580Е-01 
*696E-01 
«659Е-01 
«615Е-01 
*579E-01 
*552E-01- 
«51ЗЕ-01- 
*672E-01- 
«686Е-01 
*672E-01 
“722Е-01- 
“702Е-01 
“643Е-01- 
+56 2Е-01 
*S33E-01 
•596Е-01 
+548Е-01 
+558Е-01- 


Бр 

ЕТ/5ЕСО 
«000Е%01 
+ 433Е+02 
%197Е%03 
«395Е%03 
«551Е%03 
+616Е+03 
•556Е+03 
•195Е+03 
*142E*03 
266Е+02 
+8 OSE+03 
864E*03 
*974E*03 


.227Е+03 
175Е+03 
*140E*03 
562Е+03 
+8 43Е+03 
101Е+04 
• 103Е+04 
*911E*03 
*T754E* 03 
*384E*03 
“79 3E* 02 


103E*04 
*941E*03 
*363E«03 
*271Е+03 
164Е%03 
137Е+03 
•545Е+03 
BA2E*03 
•101Е+04 
• 103Е+04 
920Е+03 
*761E*03 
.360Е+03 
•791Е+02 


THTA 
RAD 
“000Е%01 
%507Е-01 
•251Е+00 
*583E«00 
+ 102E+01 
152Е+01 
«206Е%01 
*262Е+01 
+322E+01 
*244Е+01 
4178Е%01 
+ 109Е+01 
+360E+00 
+376E+00 
4110Е%01 
4191Е%01 
«244Е%01 
“301Е%01 
“361Е%01 
«203Е%01 
*133Е+01 
“588Е%00 
+ 166Е+00 
*906E«00 
4162Е%01 
.229Е+01 
“99 5Е+01 
%364Е%01 
4196Е%01 
+123E+01 
“873Е%00 
+28 4E+00 
+ 103E+01 
*185E*01 
%239Е%01 
*296Е+01 
*357E«01 
+206E+01 
+136Е+00 
+ 622Е+00 
*134E*00 
.877Е+00 
•159Е+01 
.226Е+01 
*347Е+01 
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ALPHA 
DEG 
«000Е%01 
4113Е%01 
%557Е%01 
+128Е+02 
*216Е+02 
+307E+02 
+38 4E+02 
.413Е+02 
*413Е+02 
*410E«02 
%347Е%02- 
+ 229E+02= 
“797Е%01- 
«833Е%01 
+233E+02 
*367E«02 
*410Е+02 
+ 413Е+02 
+413Е+02 
+ 38 0Е+02- 
«273Е%02- 
*129Е+02- 
•370Е+01 
»194Е%02 
.322Е+02 
*404E«08 
*418E*02 
•413Е+02 


• 630Е+01 
*218E«02 
“358Е%02 
«409Е%02 
*413E*02 
+41ЗЕ+02 
+38 4Е+02- 


+28 0Е+02 
+136E+02- 
.299Е+01 
+ 188E+02 
+31 7E+02 
+ 402E+02 
+ 413E+02 


FCAM 
LB 
«000Е%01 
112Е%00 
+100Е+01 
*378E*01 
+8 76Е+01 
•149Е+02 
+18 6Е+02 
*741Е+01 
*540Е+01 
“761Е%01 
773Е+01 
492Е%01 
%127Е%01 
%397Е%01 
*944E*01 
+98 3Е+01 
4855Е%01 
• 666Е+01 
4531Е%01 
*981E«01 
«988 Е%01 
*463E*01 
*413E*01 
%130Е%02 
+ 194E+02 
*119E«02 
«302Е%01 
+ 728Е+01 
911Е%01 


«341Е%01 
.906Е+01 
“109Е%02 
•101Е+02 
*623E*01 
.521Е+01 
•965Е+01 
101Е+02 
*499E*01 
+374E+01 
4127Е%02 
4192Е%02 
•131Е+02 
<570Е%01 
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TFRICT 
LB-FT 
«000Е%01 
*252Е-01 
*541E-01 
«780Е-01 
*951E-01 
•105Е+00 
%109Е%00 
*116E*00 
4123Е%00 
•127Е+00 
+132E+00 
+ 143E+00 
+ 147E+00 
*147E*00 
4144Е%00 
“102Е%00 
%108Е%00 
*117E*00 
*125Е+00 
*134E*00 
*145Е+00 
•150Е+00 
•150Е+00 
•145Е+00 
+138E+00 
+132Е+00 
4135Е%00 
%136Е%00 
%143Е%00 
*149E*00 
+151Е+00 
•151Е+00 
•147Е+00 
•105Е+00 
•110Е+00 
*118Е+00 
•126Е+00 
%134Е%00 
*145E*00 
+ 151E+00 
*151E*00 
*146E* 00 
%138Е%00 
%132Е%00 
+136E+00 


LEES 
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EPAM PROGRAM LISTING 


EFAM MATH MODEL 
VERSION OF MAY 3, 1971 


COMMON/FORCES/KI»FI^FIg,IEGUI Vs TFRI CT» MASS» RFM» РОМЕН» DDY > 

б» FN» ТРЕСТ» TH TA» ALFHA» TORQUE» TCAMs C1» C2 СЗ» CA KS» УТ, Y BOLT, 
NCYCLEs THETA, IMTEST? I FIN 

COMMON/DATA1/DATACB6) 

DIMENSION TSAVE( 250), ҒЫ5С 250)» TBELTSC250) » DDY SC 250)» THTASC 250), 
ALFHASC 250) » FCAMS( 250) › TFRICSC250)  THC 10)  REMTC 20) » ТТАВС 20) » 
ТНКАМС 3) » THPING 4) 

REAL IEGUIU^IGEAR,IFEED^I САМ» KT» MASS KS» КР» КЕ» КРІМ, МВ»МС»1.САМ 
DOUPLE PRECISION DT 


EQUI VALENCE CDATACI) ›АСАМ 2» (DATACO) »GFIELD)» 
СРАТАСЗ) »RADIUS)» (DATAC4) »КАр5РЕ), CDATACS) „СК ›, 
СРАТАС 6) ›ОМЕСА >» (DATACT) ,МС )» СОАТАСВ) „МВ ›, 


CDATAC9) sWTBELT)s (DATACIO)LCAM 9 


10 


11 


12 
13 
14 
15 
16 
17 
18 
19 
20 
21 
22 
23 
24 
25 
26 
27 
28 
29 
30 
31 
32 


33 
34 
95 
36 
37 
38 

39 
40 
4l 

42 
43 


44 
45 
46 
47 
48 
49 
50 
51 


с: 


с: 


с: 


с: 


са 


10 


п 


EGUI VALENCE, CDATACIID2IGEAR ), (DATAC12),1FEED 9» 
СОАТАС13)›1САМ 2)» СРАТАС 14) ›КР э, СПАТАС159КТ >, 
СРАТАС 16)» КЕ 9» CDATACIDKRIN 5» CDATACIBDoYINIT 2» 
СОАТАС19)›ҮРКЕ 2» (ПАТАС2О),А5ІМ > 

ECUI VALENCE (DATAC21)4THIMP >» (DATAC22) , THREL 9, 


CDATAC?3) › THTCRUD > 


СРАТАС 24) › ТНКАМС I2) 4C РАТАС 27), ТНЕТМС 1)» 


(DATAC31)4 ТНС1> 2» (ПБАТАС38),Е5ТЕР1), СРАТАС 39 > 4 FSTRF2D ^ 
€ DATACAO) › FLOCK1) 


EQUI VALENCE € DATACA1) UR э, CDATAC 42) UT ›, 
< РАТА‹ 43) › ВВ )» СБАТАС 44) CB э, (DATACAS) » DB ›, 
(DATAC 46) EB 2» CDATACAT) »RFMTC 12 > 4 DATAC 67)» TTABC ID 


FI=3+141593 FI2-2.*Flj 6=32.17 

READ INFUT DATA 

KNO=0 

АССЕЕТ МО, DTX» DX» TFINALsKPRINTs CYCLES 

KNO=KNO+1 

IFCKNO.LE.ND)CACCEFT 1, АТАСІ); GO TO 6) 

KS=1+0/€1+/KF+16/KR) 

TEQUI V=I CAM+I GEAR+I FEED. 

MASS=MB+MC 

THOS) =PL2+THC1)3 THCÓ)SFIZ*THC2)3 THC 7)=F1 2+THC 3) 
CALCULATE COEFFICIENTS FOR NORMAL FORCE ECUATION 

C1SUR+UT*( 1.+EB/5B) *2«* DP/CB*UT-2+*DE/BE*URUT 

C2-UR*UT*( 1-+EB/EB) *2 & DB/CB&UT*UR+2+* DE/EBKUT- 1.0 

C3=-(UR+UT*( |» *EB/BE) *2 & DE/CE*UT* 2 * DB/BBKÜE*UT 

C4sUR*UT*C | + +EB/BB) +2.* DE/CB«UT&UR- 2. DB/BBKUT- 1+0 
CALCULATE IMFACT FACTOR FOR ВОЦТ-САКЕТЕК IMPACT 

DYs2.*ACAMKEADIUS*CTHIMF-FI2) 

ALFHA=ATANC DY) 

REFEC= DY *KADI US**2 

VIMPAC=( (MC+I EQUI V&SINCALFHA) *COS( ALFHA) /KEFEC) - СЕЖМЕ) 
/(MC*IEGUI UKSINCALFHA) *COSCALFHA) /REFEC+MB) 

Те0.; DDX*0.0; NOUT=1 

KI=03 1№=0; Кеј) KTEST=03 LTEST=03 ТМТЕЗТ=1; KSTEDY=0 

FN=0+03 FOWER=0+03 FMUX=0+03 FOLD=0+03 IFIN*O 

Аһя0.0 3 АКТ=0.0; XOLD=0+03 XOLDT=0.03 ҮТ=0.0 
WRITE HEADING AND INITIAL CONDITIONS 

WRITEC 1430) Ts DT, TFINALsKFRINT 

WEITEC 1,60) ICAM» I GEAR» I FEED» RADI 05» MASS 


WEI TEC 1» ЗЭК, Ts Xs IKs рх, TTAEC 19/1926 РОМЕН» EE M+ Y BOLT# 126 
BEGIN INTEGRATION CYCLE. 


LFM. ЕО. 1 «О M4 ЕС. 20k «MM ЕС. 4 ORM. 2046) GO TO 11 
DAReCFOLD*POWER)*(X-XOLDT)/2«3 ARTEART*DAR; 
XOLDT=X3 FOLDstOWER 

CALL ACCEL РХ» Т»Х › DX» IN» ЕСАМ) 
ТЕСМ.ЕС.1.ОБ-М.ЕС.2.ОЯ.М.ЕС+4-ОН.М.Еб.6> GO TO 24 
KTEST#KTEST+1 

FOWER=TOROQUE* DX/550+ 

IFCX «LT+2+*CYCLES*PI 2+XOLDT) GO TO 14 

FWFSSAR/CX -XOL DT 

NouT=2 

GO TO 16 


55 


52 14 
53 

54 16 
55 

56 

57 

58 

59 

60 

61 с: 
62 20 
63 

64 

65 

66 

67 


89 28 


92 
93 30 


94. 35 


95 40 
96 50 
97 60 


98 
99 


100 


101 
108 
103 
104 


HTC-AD 71-59 


ТЕСКТЕ5Т „МЕ. KFRINT)GO TO 24 

КТЕЗТ=0; К=К+1 

WRI TEC 1» 40) К» Ts ТНЕТА1*57, 29 578 » DX › DDX» TORG UEs FOVER ^ НЕМ» Y BOLTE 1 2. 
ТЅАМЕСК) =Т; FNSCK)=FN3 TBELTSCK)=TBELT3 DDYSCK)=DDY 

ALFHASCK) =АСЕНА*57.29578; ТНТАЅСК) =ТНТА; FCAMSCK)=FCAM 
TFRICSCK) «ТЕРІ 
ТЕСТ.СТ.ТЕТМАЬ GO TO 20 

IFCX «GT+THTCRU+AND+RPM«LT+20) GO TO 20 

60 TO (24.20),/NOUT 

AT 20, TERMINATE CASE AND GO TO 5 (NEXT CASE) 
FWETSART/XOLDT 


DISPLAY" " 

DISFLAY" ^ AVERAGE TRANSIENT FOWER s", ЕШЕТ, HE" 

DISFLAY" ^ AVERAGE STEADY STATE POWER s'",EWRS»" НЕ" 

DI SFLAY RE 

DISFLAY" TIME FN TPELT DDY THTA ALPHA ", 
"FCAM ТЕРІСТ" 

DISFLAY" SEC ъв LP-FT ҒТ/5ЕСО RAD DEG ", 
"LB 


DO 22 IR=1,K 

WRITECI^SO) TSAVEC IR)» FNSCI R)» TBELTSCI E) » DDYSCI F)  THTASCI RO > 
ALFHASC IR) » FCAMSCIK)» TFHICSCIR) 

GO TO 5 

АТ 24 THRU 30 UFDATE TIME, DISFLACEMENT AND VELOCITY,» 

THEN START NEXT INTEGRATION CYCLE 

IFCTHETAL+LT+THIME „ОН. ІМТЕ5Т.ЕС.1) GO TO 25 

IMTEST=1 

18-0 

DX =VIMPAC DX 

KTEST=0 

ASSIGN 10 TO NL 

IFC DK «LT -ОМЕСА- 0 

KSTEDY=1 

XOLD=XOLDT 

IFCKSTEDY.NE.1) GO TO 28 

ASSIGN 11 TO NL 

ІҒСМ.ЕС.1.08.М.ЕС. 2.О-.М.ЕС.4.ОБ.М.ЕС.6> GO TO 28 

FOLD+FOWER) *(X=XOLD) /2. 


KSTEDY+EG+1) GO TO 26 


FOWER 
TSAINDV( Т, DT» Ms INS NC) 
X=DENUCXs рх, DT» Ms NC) 
DX = DF NVC DX > рух» рт» м, NC) 
GO TO №610, 11» 
FORMATC//,25X, 9HINITIAL CONDITIONS /4X,5H TIMEI3X,9HSTEF SIZE 
BX» 10HFINAL TIME, 6X14HFRINT INTERVAL/F10«3^2F18 + 4, 14X 13) 


FORMATC/2X » ІНІ» 3X, 20HTIME THETA ОТНЕТА ›2Х» 7HDDTHETAs 
3X» 31HTOROUE POWER RPM BOLT /6X»+12HSEC DEG 
23X233HRAD/SEC RAD/SEC2 LB-FT НЕ 9 14Xs2HINS 


132 F7 + 3» 2Р9 2, Е9 + 1, 2F9 + З» ЕВ + 1, F8 3) 

FORMATCI 3, Е7 + 3» 2Р9 + 2» F9 1» 2F9 «3s FB + 1» F8 +3) 
FOFMATCF6« 3, 7E9 +3) 

FORMATC/10X» 1BHMOMENTS OF INERTIA „9Х,ЗНСАМ 6X 4HBOLT/ 
5X, ЗНСАМЕХ SHGEARS 2X 9HFEED SYS 6X 6HRADIUS АХ 4HMASS/ 
1gX^10HFT-LB-SEC? 16X 2НЕТ 6X SHSLUGS / 

2Х,3Ғ9.7,4Х 259.5) 

END 


SUBROUTINE ACCEL( DDTHET» T» THETA» DTHETAs I Ns FCAM) 


COMMON/FORCESZKI › РІ, Fl 24 IECUI V» TFRI СТ» MASS ^ RPMs РОМЕН» DDY » 

б» FN» TBELT^ TH TA» А. РЕНА, TORGUE» ТСАМ, C1» Се, СЗ, CA» KS» Y T YBOL T^ 
МСУСЬЕ» ТНЕТА1»1МТЕ5Т,ІРІМ 

COMMON/ DATA1/DATAC8 6) 

DIMENSION THC10)^ RFMTC20), TTABC 20)» ТНВАМС 3) › THEINC 4) 

REAL IEQUI Və I GEAR» I FEED» I CAM» КТ» MASS» К, KFs KR» КРІМ» МВ» МС» LCAM 


EQUI VALENCE СРАТАС1) ›АСАМ 2, CDATAC2) »GFIELD)^ 
СРАТАСЗ) »RADIUS)s (DATACA) »RADSPR)» СРАТАС5) „СК de 
СРАТАС6) „ОМЕбА 2, (DATACT) ,МС 3» €DATACB) „МБ de 


СОАТАС9) sWTBELT)» (DATAC10)sLCAM › 


105 


106 


107 


108 
109 
110 
111 


112 
113 


сї 
Са 
Ct 
са 


с: 
С: 


а Ct 
5 с: 


116 


117 
118 
119 
120 
121 
122 
123 
124 
125 
126 
127 
128 
129 
130 
131 
132 
133 
134 


135 
136 


137 
138 
139 
140 
141 
142 
143 
144 
145 
146 
147 
148 
149 
150 
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EQUI VALENCE CDATACII)^IGEAR 2» (DATAC12)s1FEED 9, 
CDATACI3)4ICAM 2, (DATACIA) 4KP )» СОАТАС159,КТ >» 
C DATACI62 ^ KR 3» (DATAC17)s4KPIN >» €DATACI8)sYINIT 9» 
СРАТАС19),УРКЕ 2)» ¢DATAC20)sASIN 9 

EQUI VALENCE ( DATAC 21) THIMP >» (DATAC22)sTHREL 2» 


€DATACE3), THTCRU) ^. ( DATAC24) › THRAMC 122 › C РАТАС 27) › ТНРІМС 12)» 
C(DATAC31), THC1) 9» СРАТАС 38 >›Е5ТЕР!>» (DATAC39)sFSTRP2)s 
СРАТАС 409» FLOCK1) 


EQUI VALENCE, € DATACA1) ^ UR 9» (DATACARD UT >» 
€ DATACA3) „ВЕ э, (DATACAA4) › CB. )» (DATACAS) » DB. ›, 
СРАТАС 46)›ЕВ >» CDATACAT) 4 REMTCI22 5 CDATAC6D ^ TTABCID > 


DEFINITION OF BUILT-IN CONSTANTS 
*TA2869MAXIMUM SLOPE OF PARABOLIC PORTION OF CAM 
+12703=TOTAL CHANGE IN SLOFE OVER CYCLOIDAL FORTION OF САМ 
+198 4=ТОТАЬ AXIAL DISFLACEMENTCFEET) AT END OF CYCLOI DAL 
PORTION OF CAM 
*2747=TOTAL AXIAL DISFLACEMENTCFEET) AT END OF CONSTANT 
VELOCITY PORTION OF САМ 
*377210TAL AXIAL DISPLACEMENTCFEET) WHERE FIRING PIN SPRING 
COCKING BEGINS DURING REARWARD MOTION 
+0088=PRELOAD ON FIN SPRING( FEET) 
«829» TOTAL AXIAL DISPLACEMENT WHERE FIRING PIN SPRING 
CONTACT 15 MADE DURING FORWARD MOTION 
+0375=DEFLECTION OF FIRING FIN SPRING WHEN CONTACT 15 
FIRST MADE DURING FORWARD MOTION 
„04533“МАХТМОМ DEFLECTION OF FIRING FIN SPRING 
RPM=DTHETA/( PI¥2+) 60+ 
CALL TLKUPICRPMTs RPMs 20^ KI» AY Y 05 0» 1, TTAB» TORQUE» ER) 
10-21; АРАСТ=1.0; TORQUE=TORQUE/192+3 FLOCK=0+0 
ТНТА=ТНЕТА; THETAISTHETAS ҮРІМе0.0) FSTRIP=0+03 MASS=MC+MB 
IFCTHETAsLT+2+*F12) GO TO 2 
МСУ CLE=THETA/( 2.*Р1 2) 
THTA«THETA-2.*NCYCLE«PI? 
ТНЕТА1=ТНТА 
ІЕСТНЕТА1.6Т.ТНЕЕІ..АМО»ТНЕТА1.1.Т.ТНІМР) МАЅ5=МС 
ТЕСТНЕТА! «LT+THRAM( 1) -OR«THETAL+GT+THRAMC2)) GO TO 4 
FSTRIPAFSTRFI 
ТЕСТМТЕЗТ. ЕС. 1) 1№=0 
IMTEST=0 
IFC THETA] • СТ. THRAM( 2) AND. THETA +L T+ THRAM( 3) ) FSTRI P= FSTRP2 
IFCTHTA+GT+TH(3)) GO TO 6 
COF=1+ 
YBOL T=ACAM*( THTA*RADI US) **2 
GO TO 14 
IFCTHTA+GT+FI2) GO ТО 8 
ТНТА=Р12-ТНТА 
С0Е=-1.; 1082 
YBOLTSLCAM-ACAM*(THTA*RADIUS)*«2 
IFCTHTA«LE« THCI)) GO TO 28 
GO TO 15 
IFCTHTA «GT. THC72). GO TO 10 
THTA®THTA=PI 2 


1023 
YBOLTSLCAM-ACAM*(CTHTAKRADIUS)**2 
GO TO 18 


THTAs2.*PI2-THTA 

COFs-1.; 16=4 

YBOLT=ACAM#( THTA®RADI US) 4*2 

GO TO 7 

IFCLQsEG+1 «AND.THTA.LE. THC1)) GO TO 23 
IFCTHTA+GT+TH(2)) GO TO 16 

GAM=( THTA~THC 1) )*PI7( 2" ASIN? 
АҒАСТ=С05‹ GAM) 

DY=+74226++12703*S1NC GAM) 
1Ғ<10.Е0.1.08.16.Е6.3) GO TO 24 

GO TO 29 

СОЕ=1.0 

YBOLT2.138 4+( THTA- THC 2) ) RADI USKTANC AL FH A) 
GO TO 22 

IFCIG.NE.3 „ОВ. ТНТА.бЕ.ТНС1)> GO TO 20 
СОҒ=-1.0 

60 то 23 
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164 20 IFCIG+NE+3 «ОК. THTA«GT«THC2)) GO TO 21 


165 GO TO 15 

166 21 COFs-1. 

167 YBOLT= +2747~¢ THTA~TH( 2) ) RADI US&TANCALFHAD 

168 Ct STATEMENTS 22 THRU 27+2 ARE ACCELERATION EGUATIONS 
169 22 ALFHA-.72 

170 DDY «RADI US*DDTHET* TANC ALPHA) 

171 GO TO 25 


172 23 DYs2.«ACAM*RADI USKTHTA 
173 24 ALPHAS=ATANC DY) 


174 Юру =2.*АСАМ*( RADI US** 2) ж DIHETAW 2) XAF ACT+RADI USKDDTHET* DY 

175 25 ТЕСТНЕТА. GT THTCRU). Y TeY PRE*Y BOLT-YINIT 

176 ІЕСТНЕТА1.1,Е.2.4РІ2-ТНРІМС2).ОБ «IFIN«NE«1) GO TO 42 

177 IN=0 $ ІРІМ=0 

178 42 ІЕСТНЕТА1.1.7.2.ж-12-ТНРІМС 3) «OR! THETAL бТ.2.*Р12-ТНРЕМС 222 
60 TO 26 

179 ІРІМ=1 

180 YFIN@ACAM*( TH TA®RADIUS) **2-CLCAM- • 377- • 0088» 

181 FLOCKs1«5*FLOCK1 

182 26 BFMASS*DDY+KS*YT+KFIN#*Y PIN+FSTRI F+FLOCK 

183 IFCBF«LT+0+) 60 TO 31 

184 27 FN=-BF/(C1*SINC ALFHA) +C2*COSC ALFHA)) 

185 FCAM=FN#( SINC ALPHA) +UR*COS¢ ALPHA) У 

186 DDY =DDY*COF 

187 GO TO 40 

188 31 BF=-BF 

189 GO TO 36 


190 32 BF=-BF 

191 GO TO 27 

192 C: STATEMENTS 28 THRU 36+3 ARE DECELERATION EQUATIONS 
193 28 DYz2.*ACAM*RADIUS*THTA 

194 29 ALFHA=ATANC DY) 


195 DDY = 2.жАСАМЖ RADI 05жж 2) «C DIHETA®*2) AFACT-RADI USKDDTHET& DY 

196 IFC THETA] «LT. THEINCI) „ОВ ТНЕТА1 •67.ТН РИМ 2) ) 60 ТО 30 

197 YFINs.0375*LCAM- „ 329 ~ACAM*( THTAKRADIUS)**2 

198 FLOCK=FLOCK13 UCOEFs1«0; GO TO 35 

199 30 ТЕСТНЕТА! «-Т.ТНРЕМ 2) „ОК. THETA] • GT THEINC 37) 60 TO 35 

200 ҮРІМа.04533 

201 FLOCK=FLOCK1; UCOEFsO.3 

202 35 BF=MASS*DDY -FSTHI F-FLOCK-UCOEF*KPIN&YPIN 

203 IFCBF.LT.O.) GO TO 32 

804 36 FN=-BF/(C3*SINC ALPHA) +C4*COS( ALPHA) > 

205 FMUX = UR*FN#COS¢ ALPHA) 

206 FCAM=FMUX -ЕМ+51 С ALPHA) 

207 DDY = рру *COF 

208 40 CONTINUE 

209 TPRELTsWTBELT&RADSPRAC1«B*GFI ELDEDDTHET*RADSPR/CI2«*G)2/12« 

210 TCAMSFCAM*RADIUS 

211 TFRICT=KT*DTHETA 

212 DDTHET=( TORQUE- TCAM- TFRI CT- TBEL T) /I EQUIV 

213 RETURN 

214 END 

215 BLOCK DATA 

216 COMMON/ DATA1/DATA?( 10)» РАТАЗС 30) » DATAAC 6) › DATASC 40) 

217 DATA ПАТА2/1.6609»1.0» «11717» «0625» „75, 185.6667, «017811» 
0118744» 6.4» «413087 

218 DATA ПАТАЗ/.0000277» „00000165, «00087» 740»».001,162852.» 


528.» «0883, «045545585» 7.56564» 5.0, 100», 3.264» 3» 38 6» 
4.433) 4.3633, 4 459 3» 4.9742, Os» 1.9199, 2.478 4» 3-8 048 › 4- 3633» 
0%,0.)0.,20.0,6.8,2./ 

219 DATA БАТА4/.063» „195, 2.54, 1.72» «99, .58/ 

220 DATA ОАТА5/0.0»100.,200.»300.,400.»500.,600»» 700.>800.,975.» 
1037.5,1112.5,1175.»1237.5» 1287.5» 13504» 1400 +» 1460», 160002170008 
1120.»1110»»1100.»1070.» 1030.,985.,930.2865.,785.,640.»560»,480.» 
400+» 320.»240.,160.,80.,0.,-320.»-600./ 

221 END 


38 
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222 С:ОРУМ SUB. TO INTEGRATE BY ADAMS METHOD 
223 FUNCTION ALNDV( Ts рт, M> ININ) 

224 DOUBLE PRECISION DT 

225 Nel 

226 IFCIN „МЕ. 0260 TO 2 

227 1 м=о 

228 INs1 

299 2 МаМФі 

230 IFCM=7) 3» 306 

231 3 GO T0(6,5555 6/5565 55, M 

232 5 AINDV=T 

233 RETURN 

234 6 AINDV=T+DT 

235 RETURN 

236 END 

237 FUNCTION DPNVCY S» DYS» DT, MN) 

238 DIMENSION YC20), DYC4^ 20)» DY11¢20) 
239 DOUBLE PRECISION DT,Y,DY,DYlloYY 
240 YYsYS 

241 DY Da DY S 

242 ТЕСМ-728»8»9 


243 8 GO то с10,11,12,13, 1415» 16),М 
244 10 DYC2sN)=DYD 


245 YOND=YY 

246 XYYeY CN) +DT* DY D 

247 17 DPNV=XYY 

248 N=N+1 

249 RETURN 

250 11 DY11(N)=DYD 

251 хумус *CDT/2 4) *C DY 1168) руса, ND 

252 60 TO 17 

253 12 KYY=YCN)+( DT/6+) «CO e* DY D+ DY 1161) *3«4XDY CONDO 
254 GO TO 17 

255 13 Y(M=YY 

256 0ҮС 3, № =рүр 

257 ХҮҮзҮСМ)жСВТ/2.3жС 34 DY C34 ND -ру Се» NDO 

258 60 TO 17 

259 14 XYYsY(ND)*CDT/ 19.) «C5  DY 0+8 DL C 3» ND -DY CONDO 
260 GO то 17 

261 15 YCND-YY 

262 рус 4,N)=DYD 

263 XYYsYCN) *CDT/12.2 «C 236%DY D- 16«* DE C34 №+5.ж*1 (2, №) 
264 GO TO 17 


265 16 XYY=YCN)+( DT/ 24+) *C9 «DY DE 19 DYCASND -5«* DYC 34 ND «үс 2,899 
266 GO TO 17 
267 9 DO 18 I*1,3 


268 IIsI*1 

269 18 DYCI,NDeDYCII,N) 

270 DYC4,N)=DYD 

271 XYYSYY+CDT/ 246) CSS Dr C4» №) -59 DY C3 № +37•жру се» ND -9 «жс, МӘ) 
212 GO TO 17 

273 END 

274 SUBROUTINE TLKUP! (AK sXsLXsKe AY »Y»LY» Jo NO» АС» ANS» ER) 

275 DIMENSION АХС 80)» АСС 20) 

276 300 ER = 0. 

277 IFCK=1) 310» 310, 311 


278 311 ТЕСАКСК)у-Х) 316, 313» 313 

279 313 ТЕСАКСК-1)-Х) 390, 317, 308 

280 317 ТРСК-2)312, 390» 312 

281 308 IF (K-2) 312» 307» 312 

282 307 МРІТЕС1,998) X» АХСКО 

283 998 FORMAT ( ІН 48HEXTRAPOLATION LOW END, 2 D TABLE, VARIABLE 1= 
Е15.6» АКЛАНТАВЊЕ VALUE = Е15.6 > 


284 60 TO 329 
285 310 K= 2 
286 GO TO 311 


287 316 IF(K-LX)3155309,310 
288 312 K = K-1 


289 60 TO 313 
290 315 K = ке 
291 60 TO 311 


39 


292 
293 


294 
295 
296 
297 
298 
299 
300 


309 
999 


329 
390 


395 
403 


HTC-AD 71-59 


WRITEC 1, 999) X, AXCK) 

FORMAT € 1H 48HEXTRAFOLATION HIGH END, 2 D TABLE, VARIABLE 1= 
Е15.6» АХ1АНТАВЬЕ VALUE = Е15.6) 

ER = 140 

ХМ = (Х-АХСК- DD) /САХСК) - AX CK 12 

DO, 395 KOUNT=1»NO 

M = LX ж CKOUNT - 10 + K 

ANS = ACCM-1) + XM*CACCM) = ACCM=1)) 

RETURN 

END 


40 


